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ABSTRACT 

While stars are widely discussed as the source of the high energy photons which reionized the universe, 
an additional source of ionizing photons that must also contribute to reionization in this scenario is 
the supernovae (SNe) which mark the end of the life of massive stars. Here we estimate the relative 
contributions of SNe and stars to reionization. While the rate at which ionizing photons are produced 
in SN shocks is well below that at which they are produced by stars, the harder spectra of radiation 
emitted from SNe leads to an enhanced escape fraction of SN generated photons relative to that of stellar 
photons. In particular, along a given line of sight out of a galaxy, we find that for neutral hydrogen 
column densities iVn > 10^* cm~^ the contribution to reionization from SNe is greater than that from 
stars. Drawing on the results of simulations presented in the literature, we find that the overall (line 
of sight-averaged) SNe shock-generated ionizing photon escape fraction is larger than the stellar photon 
escape fraction by a factor of ~ 4 to ~ 7, depending on the metallicity of the stellar population. Overall, 
our results suggest that the effect of SNe is an enhancement of up to ^ 10 percent in the fraction of 
hydrogen reionized by stellar sources. We briefly discuss the implications of our results for the population 
of galaxies responsible for reionization. 

Subject headings: supernovae - reionization 



1. INTRODUCTION 

The reionization of the universe was one of the 
paramount milestones in the transition from the cosmic 
Dark Ages, the period before the formation of the first 
stars, to the present day (see e.g. Barkana & Loeb 2001; 
Ciardi & Ferrara 2005). As it is the first generations of 
luminous objects that were the sources of the high en- 
ergy radiation which reionized the intergalactic medium 
(IGM), theoretical studies of how reionization took place 
(e.g. Gnedin & Ostriker 1997; Miralda-Escude et al. 2000; 
Iliev et al. 2006; McQuinn et al. 2007; Shin et al. 2008) 
in conjunction with observations of the radiation emitted 
during the epoch of reionization (e.g. Fan et al. 2006; 
Robertson et al. 2010; Bowman & Rogers 2010) stand to 
teach us a great deal about the nature of the first stars 
and galaxies (see e.g. Bromm & Yoshida 2011). 

One of the outstanding questions pertaining to reion- 
ization is what were the sources of the ionizing radiation 
which effected it. It appears unlikely that the progenitors 
of the galaxies observed to date at z > 7 are responsible for 
reionizing the universe, as the escape fraction of ionizing 
photons from these galaxies would have to be very large or, 
otherwise, a low-metallicity stellar population or one with 
a relatively top-heavy mass function would be required 
(see e.g. Ouchi et al. 2009; Bunker et al. 2010; Wilkins 
et al. 2011; but see also Wyithe et al. 2009; Bouwens et 
al. 2011). However, there are a number of possible sources 
which likely contributed to reionization at higher redshifts, 
where observations of the cosmic microwave background 
(CMB) suggest a large fraction of the IGM was reionized 
(Komatsu et al. 2009). Among these are stellar sources, 
including possibly massive Population (Pop) HI stars (e.g. 
Ciardi et al. 2003; Yoshida et al. 2004); quasars pow- 



ered by accretion onto black holes (e.g. Ricotti & Ostiker 
2004; Alvarez & Abel 2007; Volonteri & Gnedin 2009); 
early populations of high-mass X-ray binaries (e.g. Power 
et al. 2009; Mirabel et al. 2011); high-velocity structure 
formation shocks (Miniati et al. 2004; Dopita et al. 2011); 
and dark matter decay and annihilation (e.g. Schleicher 
et al. 2008). 

One additional source of hydrogen ionizations which in- 
evitably accompanies stellar sources are the supernovae 
(SNe) which mark the end of the life of massive stars. 
The impact that SNe have on reionization is at least two- 
fold: not only are ionizing photons generated in SNe ex- 
plosions (e.g. Chevalier 1974; ShuU & Silk 1979), but 
the coUisional ionization and evacuation of gas they cause 
lower the optical depth to ionizing radiation, thereby al- 
lowing for a higher escape fraction of ionizing photons into 
the IGM (e.g. Clarke & Oey 2002; Yajima et al. 2009; 
Paardekooper et al. 2011; see also Tegmark et al. 1993). 
At the high redshifts {z ^ 20) at which the first stars 
formed, due to the high energy density of CMB photons 
(which is oc {\+z)'^) X-ray photons produced via Compton 
scattering of the CMB inside SNe remnants are likely to 
have contributed to the reionizization of the IGM, as de- 
scribed by Oh (2001). However, over the full redshift range 
over which reionization takes place (i.e. z > 6) SNe gen- 
erate ionizing photons predominantly through other pro- 
cesses, such as bremsstrahlung and resonance and recombi- 
nation line emission (e.g. Cox 1972; Chevalier 1974; Shull 
& McKee 1979); indeed, these processes can also be dom- 
inant even for some Pop HI SNe in minihaloes at z ~ 20 
(see Kitayama & Yoshida 2005; Greif et al. 2007; Whalen 
et al. 2008). While, in any case, the number of ionizing 
photons produced in SNe shocks is in general well below 
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the number produced during the pre-SN phases of stellar 
evolution, because they are on average more energetic than 
the photons emitted from stars, the cross section for their 
absorption by atomic hydrogen and helium is smaller; this, 
in turn, allows for a higher fraction of SN shock-gcncratcd 
photons to escape into the IGM and contribute to reion- 
ization. 

In the present work, we calculate the magnitude of this 
effect and derive an estimate of the contribution, relative 
to stellar sources, that ionizing photons generated by SNe 
make to cosmic reionization. In the next Section, we es- 
timate the number of ionizing photons generated in SNe 
shocks using recently published high velocity shock mod- 
els. In Section 3, we calculate the average cross section 
for their absorption by neutral hydrogen and helium, and 
we use this to estimate the escape fraction of SNe shock- 
generated ionizing photons as a function of the escape 
fractions of stellar photons found in the detailed numeri- 
cal simulations of Gnedin et al. (2008) and Razoumov & 
Sommer-Larsen (2010). In turn, we draw on these results 
to estimate the relative contributions of stars and SNe to 
reionization in Section 4. Though of less overall impor- 
tance, in Section 5 we furthermore provide an upper limit 
for the relative contribution of cosmic rays generated in 
SNe shocks to reionization. In Section 6 we briefly dis- 
cuss the implications of our findings for the population of 
galaxies responsible for reionization. Finally, we discuss 
our results and give our concluding remarks in Section 7. 

2. THE RATE OF HYDROGEN IONIZATION 

We begin by estimating the rates at which hydrogen 

atoms are ionized due to both SNe and stars, as a func- 
tion of the star formation rate. In the following Sections, 
we calculate the average cross section for the absorption 
of these photons by neutral hydrogen and helium, and we 
use these quantities to compare the relative contributions 
of SN-generated photons and stellar photons to the reion- 
ization of the IGM. 

2.1. Ionization Rate Due to Supernovae 

To estimate the rate at which hydrogen atoms are ion- 
ized due to the ionizing radiation emitted from SN rem- 
nants, we first estimate the total energy emitted in ioniz- 
ing photons from SN shocks during the early Sedov- Taylor 
phase of their evolution. We then estimate the average en- 
ergy of the SN shock- generated ionizing photons, and use 
this to estimate the rate at which atomic hydrogen is ion- 
ized, as a function of the star formation rate. 

2.1.1. The Energy Emitted in Ionizing Radiation 

In order to estimate the rate at which ionizing photons 
are produced in SN shocks we make use of the MAPPINGS 
III high velocity shock models presented in Allen et al. 
(2008). We use an approximation of the convenient for- 
mula (their equation 5) provided by these authors to cal- 
culate the fiux Fion of hydrogen-ionizing radiation, mostly 
composed of thermal bremsstrahlung and recombination 
and resonance lines, that is emitted from shocks with ve- 
locities 10^ km < Vsh. < 10'^ km s~^ passing through 
a medium with a particle number density 0.01 cm~'^ < n 
< 10^ cm-^: 



cm~ 



erg cm ^ s ^ 



Next, we make the assumption, justified below in Section 
2.1.2, that ionizing photons arc produced by SNe shocks 
while in the Sedov- Taylor phase of their expansion (see 
e.g. Truelove & McKee 1999). Using the Sedov- Taylor 
solution, we thus find the following expression for the the 
distance rsh that the shock has traveled since the initial 
SN blast, as a function of its velocity Vsh- 



rsh = 35 



\ 10^1 erg 
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(2) 

Using equations (1) and (2), we integrate over the whole 
of the spherical blast wave and over time, to find the total 
energy Eion emitted in ionizing radiation : 
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where in the second part of the equation we have expressed 
the integral in terms of Vsh using equation (2). The lower 
limit of integration is set by the approximate minimum 
shock velocity (~ 200 km s~^) for which ionizing radiation 
may propagate beyond the shock, assuming a largely neu- 
tral medium ahead of the shock (sec e.g. ShuU & McKce 
1979; Allen et al. 2008); for lower velocities, the ionizing 
photons produced in the shock are all absorbed by neutral 
gas as it is swept up by the shock. If the ambient medium 
is instead largely ionized, then ionizing radiation produced 
in somewhat lower velocity shocks will propagate beyond 
the shock as well; thus, our assumption of a neutral am- 
bient medium is a conservative one (see e.g. Shull & Silk 
1979). The upper limit of integration in equation (3) is set 
by the maximum velocity for which the shock models of 
Allen et al. (2008) are given; however, for this reason, too, 
our estimate of the energy emitted in ionizing radiation is 
conservative. 

Also, interestingly, there is no dependence on the density 
n of the ambient gas in equation (3). Rather, it is simply 
that ^ 40 percent of the energy of the SN is emitted as ion- 
izing radiation, which is roughly consistent with previous 
estimates (see e.g. Mansfield & Salpeter 1974; Chevalier 
1977), as well as being consistent with the amount of en- 
ergy that is radiated during the Sedov- Taylor phase (see 
e.g. Draine & Woods 1991; Shu 1992), the extent of which 
we shall turn to discuss in the next Section. 

Before moving on, we must note that the shock models 
presented in Allen et al. (2008) arc, strictly speaking, valid 
in the steady-state approximation, which may be a poor 
description of the early stages of the Sedov- Taylor phase. 
Nonetheless, we find agreement with previous works in 
part because the majority of both the total energy radi- 
ated in ionizing photons and the total number of ionizing 
photons (see Section 2.1.3) are emitted at the late stages 
of the Sedov- Taylor phase, when the steady state approx- 
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imation is more robustlf] Furthermore, we note that even 
the gas shock-heated to very high temperatures > 10^ K 
at early times will eventually cool radiatively, albeit on 
timescales of up to ~ 10^ yr, considerably longer than the 
extent of the Sedov- Taylor phase (see e.g. Cioffi et al. 
1988). Therefore, while it may take upwards of a million 
years, the post-shock gas will cool and the time-integrated 
spectrum of the radiation emitted from the SNe remnants 
should converge, at least approximately, to the result that 
we have found under the assumption of steady-state (ra- 
diative) shocks. As this timescale is still shorter than the 
~ 10^ yr timescale on which most massive stars emit ion- 
izing photons and explode as SNe, our use of steady-state 
shock models should not greatly impact our results in the 
end. 

2.1.2. The Extent of the Sedov- Taylor Phase 

The radius rj-ad at which a SN shock transitions from 
the Sedov- Taylor phase to the radiative phase depends on 
the rate at which the shocked gas cools (see e.g. ShuU & 
Silk 1979; Blondin et al. 1998). We shah draw on the 
results of Draine & Woods (1991), who model the cooling 
of SN remnants. These authors find the following fitting 
formula for Tj-ad , for the case of a gas with solar metallicity 
(see their tables 1 & 2): 



Trad ^ 32 



VlO^^erg 



cm 



-0.45 



pc 



(4) 



For gas with a metallicity different than solar, the time 
irad at which the transition to the radiative phase occurs 
can be estimated as the cooling time of the shocked gas, 
and is thus inversely proportional to the cooling rate A of 
the gas. In turn, this affects the radius of the transition, as 
in the energy conserving phase riad oc trad^^^ oc A~^/^ (see 
e.g. Blondin et al. 1998). To express Trad as a function 
of metallicity, we can thus follow Blondin et al. (1998), 
and estimate the cooling function of the gas in the energy 
conserving phase, as a function of metallicity Z , by fitting 
to the cooling functions presented in Sutherland & Dopita 
(1993). We thus find an approximate fit to the cooling 
function as 



A ^ 10^22 ( ^ 



Z 



0.44 



erg cm'^ s ^ 



(5) 



(io^kJ 

where T is the temperature of the gas. This fit is valid 
over the temperature range ^ 10^ K < T < 10^ K, in 
which the transition to the radiative phase takes place for 
n < 10^ cm~'^ (e.g. Draine & Woods 1991), and over the 
metallicity range 0.05 Zq ^ Z < Zq, which is the same 
as that over which we model the stellar ionizing photon 
output below in Section 2.2. Using the scaling of A with 
metallicity, we follow equations (4.15) and (4.16) of Draine 
& Woods (1991) to find the following expression for ri-ad 
as a function of metallicity (for n < 10^ cm~'^): 
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Thus, there is a very weak dependence of the transition ra- 
dius on metallicity, and a somewhat stronger dependence 
on the density of the gas into which the SN explodes. Fol- 
lowing equation (2), in the Sedov- Taylor phase the dis- 
tance to which the shock has propagated, as a function 
of shock velocity Vs\i, is 



Tsh = 22 



SN 



\ 10^1 erg 



Ucm-3/ V200kms-i/ 



pc . 

(7) 

Thus, the radius at which the velocity drops to 200 km s ^ , 
the minimum shock velocity for which we assume ionizing 
photons can escape the shock front (see Section 2.1.1), is in 
general smaller than the radius at which the Sedov- Taylor 
phase ends; therefore, our assumption in the previous Sec- 
tion that the Sedov- Taylor phase lasts until Wgh = 200 km 
s~^ is sound in most cases. We note, however, that for 
-EsN ^ 8 X 10^^ erg and/or n > 40 cm"^, r^h when v^h — 
200 km s~^ will exceed r^ad and our simple treatment of 
the dynamical evolution of the SN remnant will have to be 
modified to account for ionizing photon production after 
the Sedov- Taylor phase. 

We would like to emphasize that, while the Sedov- Taylor 
solution does well to describe the dynamical evolution of a 
SN remnant at early times under the assumption of con- 
servation of mechanical energy, it is in fact the case that a 
portion of this energy is radiated away during this phase. 
Indeed, the Sedov- Taylor phase ends only when a sub- 
stantial portion of the energy of the SN has been radiated 
away, with the result that the conservation of mechanical 
energy ceases to be an appropriate approximation for de- 
scribing the dynamical evolution of the remnant, as has 
been found in myriad studies (see e.g. Chevalier 1974; 
Draine & Woods 1991; Blondin et al. 1998; Truelove & 
McKee 1999; Shu 2002). As noted in Section 2.1.1, the 
amount of energy that we have found to be lost to ionizing 
radiation, especially in the late stages of the Sedov- Taylor 
phase, is in good agreement with these previous works. 

2.1.3. The Number of Ionizations 

We estimate the total number of ionizing photons emit- 
ted per SN using the tabulated ionization parameters given 
in Table 3 of Allen et al. (2008). Carrying out an integra- 
tion similar to that shown above, approximated as a sum 
over shock velocities of Wgh — 200 to 10'^ km s~^, we es- 
timate the total number TVion of ionizing photons emitted 
per SN to be 



-EsN 



(8) 



4.4 X 10^" , 

lO'^^ erg 

Dividing -Eion by N{o-a we find the average energy of ion- 
izing photons to be < E'.y > = 51 eV, which implies that 
the majority of the energy of the ionizing radiation is car- 
ried by photons with energies > 40 eV. As shown by ShuU 
& van Steenberg (1985), ionizing photons with such high 
energies may ionize more than one hydrogen atom per pho- 
toionization, due to secondary ionizations caused by col- 
lisions with photoelectrons. In particular, these authors 
show that the fraction of the energy of such ionizing pho- 
tons that go into ionizing hydrogen is /ion ~ 0.3. We thus 



^ The reason that most energy is radiated at the late stages is two-fold: firstly, due to the remnant's slowing expansion with time, the remnant 
spends more time in the late stages; secondly, the emitting area is much larger at the late stages than at early times. 
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estimate the total number of hydrogen atoms that can be 
ionized by the radiation emitted from a single SN by tak- 
ing it that 30 percent of the energy in ionizing radiation 
-Eion goes into ionizating hydrogen atoms, each of which 
requires Eji = 13.6 eV. From this, we find that the to- 
tal number of ionizations per SN is /ion-EionZ-En — 5.2 x 
1060(i;sN/105^ erg). 

In turn, we can express the rate QsN at which hydrogen 
ionizations take place for a given star formation rate 
and a stellar initial mass function (IMF) for which A^sn 
supernovae are produced per unit mass in stars, as 

(lO-2M-i) (tO^) • 

We have normalized to A^sn = 10~^ Mg^, which is ap- 
proximately the value for a standard Salpetcr-likc IMF 
(e.g. Lcithcrcr ct al. 1999). Wc note that equation (9) 
takes into account only type II SNc produced by stars 
with masses > 9 Mq, and not type la SNe; this is reason- 
able, however, because the number of type la SNe, which 
typically occur with a significant delay following a star 
formation episode, is in general much smaller than that of 
type II SNe (see e.g. Nagashima et al. 2005). 

2.2. Ionization Rate Due to Stars 

Next, we estimate the rate Q* at which ionizing pho- 
tons are produced by stars, as a function of M^. , assuming 
a given number of ionizing photons 77 emitted per baryon in 
stars. Normalizing rjtoa value appropriate for Pop II stars 
formed with a standard Salpeter-like IMF (e.g. Haiman 
2009), we have 

The quantity 77 is metallicity and IMF dependent, and 
can, for instance, be upwards of 10 times higher than the 
fiducial value we have chosen here for a population of ex- 
tremely low- or zero-mctallicity stars (e.g. Bromm ct al. 
2001; Schacrcr 2002; Wise & Cen 2009). In practice, wc 
shall use the values of given for a range of metallici- 
ties in the Starburst99 database (Leitherer et al. 1999); 
specifically, we shall use the models from these authors for 
stellar populations formed with a Salpctcr IMF, with mini- 
mum and maximum initial stellar masses of 1 and 100 Mq. 
These choices will allow to facilitate a direct comparison 
with radiative transfer simulations from which the escape 
fractions of ionizing photons are calculated, as described 
in the next Section. 

Comparing Qsn to Q*! we see that the rate at which 
hydrogen atoms are ionized due to radiation from stars is 
in general more than an order of magnitude higher than 
that due to the emission of ionizing radiation from SN 
shocks. Nonetheless, the latter can contribute significantly 
to reionization due to a higher fraction of these photons 
escaping from galaxies, which we now turn to show. 



3. THE ESCAPE FRACTION OF SNE SHOCK-GENERATED 
IONIZING PHOTONS 

In this Section wc first calculate the average photoion- 
ization cross section for SNe shock-generated photons, and 
we then apply this to attain an estimate of the escape frac- 
tion of such photons as compared to that of stellar ionizing 
photons. 

3.1. The Average Cross Section for Photoionization 

In order to estimate the average cross section for pho- 
toionization of SN shock-generated photons, wc make use 
of the high velocity shock spectra presented in Allen et al. 
(2008). As shown in Section 2.1.3, the average energy of 
ionizing photons from SNe shocks is high enough that we 
need to account for multiple hydrogen ionizations per pho- 
ton, as wc have done in calculating Qsn- As the number 
of ionizations is thus dependent on the total energy in ion- 
izing radiation, it is the fraction of the energy in ionizing 
photons which escapes the host galaxy which determines 
the impact of SN shock-generated radiation on the reion- 
ization of the IGM. Therefore, we estimate the cross sec- 
tion CTSN for photoionization as the photon energy-averaged 
cross section: 

/i^eev i'^ME^) + ^HeUi;^)] F{E,)dE^ 

where the spectrum is described by F{Ej), the flux per 
unit photon energy. We use the photoionization absorp- 
tion cross sections for neutral hydrogen and helium pro- 
vided by Osterbrock & Ferland (2006). In particular, for 
photon energies Ey > 13.6 eV we take the cross section 
for absorption by neutral hydrogen to be cthi = 6.3 x 
10~^^ cm^ {Ej/13.6 eV)~^, whereas for photon energies 
above 24.6 eV we compute the total cross section for ab- 
sorption as the sum of the neutral hydrogen and neutral 
helium cross sections for absorption, the latter taken to 
be fJHei = 7.8 X 10-18 cm2 {E-y/2A.6 eV)-^. We have 
assumed a neutral helium to hydrogen ratio of 0.1, such 
that our results are obtained as the average cross section 
to absorption per neutral hydrogen atom, which is a con- 
venient form for the calculation of the escape fraction of 
SN shock-generated photons which follows. We note that 
this assumption implies a negligible amount of absorption 
of photons by helium which is already once-ionized (i.e. 
in the form of He ll). Given the relatively soft spectra of 
the stellar sources, the fraction of helium in the form of 
He II is likely to be low; given furthermore the low pho- 
toionization cross section of He 11 compared to that of He I 
(e.g. Osterbrock & Ferland 2006), this is a conservative 
approximation. 

The average cross section for photoionization that we 
obtain following this procedure is ctsn ^ 0.7 x 10-^^ cm^, 
for the solar metallicity shock spectra presented in Allen et 
al. (2008). As these authors note, the strength of the ion- 
izing field, being largely generated by brehmsstrahlung, is 
only weakly dependent on the metallicity of the gas; there- 
fore, we shall use this single value for our calculations for 
simplicity, although we will express our results in terms of 
crsN where possible, such that it is clear how our results 
will change for different average cross sections. 
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Next, in order to estimate the escape fraction of SN 
shock-generated photons relative to that of stellar pho- 
tons, we also need to calculate the average cross section 
<7.^, for photoionization of stellar ionizing photons. For this 
we adopt the model spectra presented in Leitherer et al. 
(1999). In particular, we use the spectra for continuous 
star formation at metallicities oi Z — 0.05, 0.2, 0.4, and 1 
Zq, and a Salpeter IMF over a mass range from 1 to 100 

Mq. 

We carry out the same exercise as described above for 
the calculation of ctsn, but we now instead calculate the 
photon numfeer-averaged cross section for photoionization. 
This is appropriate in this case, because for the softer spec- 
tra of ionizing radiation emitted from stars a single ioniz- 
ing photon will on average lead to the ionization of a single 
hydrogen atom. We thus calculate the cross section as 



given line of sight out of the galaxy does not change sig- 
nificantly between the time when a population of stars is 
shining and when they explode as SN (i.e. within ^ 3 x 
10 yr; e.g. Leitherer et al. 1999), the column density will 
be the roughly same for both cases. Assuming that it is in 
fact the same, we have for the ratio of the escape fractions 
of SNe and stellar ionizing photons 



/< 



csc,SN 
Jcsc , * 



(14) 



where = N^icr* and tsn — ^hctsn- In turn, using once 



again f^s 
photons: 



we solve for the escape fraction of SN 



fcsc, 



SN 



— /esc,* — /( 



0.27 

CSC,* ' 



(15) 



lu.GeV ['^Hl(-B^) + (THcl(-E^)] ^^^dE, 



r 

J13 



(12) 



where for the final expression we have used the values of 
(Tsn and a, found in Section 3.1 for the case of star forma- 
tion at a constant rate with a Salpeter mass function from 
1 - 100 Mq and a metallicity of Z = 0.05 Zq, as presented 
in Leitherer et al. (1999). 

We can see immediately from equation (15) that, espe- 
cially for low escape fractions of stellar photons along a 
given line of sight, the escape fraction of SNe photons can 
be much larger than that of stellar photons. For example, 
using the last expression in equation (15), for /osc,* = 0.1 
we obtain /esc,SN = 0.54, and for /esc* = 0.01 we obtain 
/csc,SN = 0.29. Taking the ratio of the number of ioniza- 
tions caused by SNe and stars found in Section 2 (Qsn/Q* 
^ 0.02), we see that more reionizations due to SNe than 
to stellar photons will take place in the IGM (i.e. that 
/esc,* Q* ^ /osc,SN Qsn) when the escape fraction of stellar 
photons is /esc,* ^ 0.005; this corresponds to a hydrogen 
column density along the line of sight out of a galaxy of 
A'h — 10^^ cm~^. For column densities greater than this, 
the contribution to reionization from SNe is greater than 
that from stars. 

While in deriving equation (15) we have assumed both 
a constant star formation rate and no evolution in the col- 
umn density along a given line of sight over a timescale of 
< 3 X 10^ yr, roughly the time between the emission of the 
bulk of the ionizing photons from stars and the emission 
from SNe remnants, in general these assumptions may lead 
us to underestimate the escape fraction of SNe-generated 
photons. The reason for this can be seen by considering 
a single stellar population formed instantaneously. In this 
case, the high energy photons emitted by stars will act to 
heat and ionize the gas along a line of sight out of the 
galaxy, and the effect of this stellar feedback will be to 
increase the escape fraction of ionizing photons with time 
(see e.g. Johnson et al. 2009; Wise & Cen 2009). Because 
this feedback occurs before the ionizing photons from SNe 
are emitted, it is likely that in general the column density 
along a line of sight out of a galaxy from a given clus- 
ter of stars is lower when SNe occur than when the bulk 
of the ionizing photons from stars are emitted. There- 
fore, it is likely the case that the escape fractions of SN- 

* In the appendix, we show expHcitly the approximation we have made in order to arrive at this formula. The reason that we do not use the 
exact formula for the escape fraction, which is the integral over E-^ of e~'^HO'(^^T) weighted by the source spectrum, is that we must extract 
the probability distribution function of the hydrogen column density dP/dNn from the escape fractions reported in Gnedin et al. (2008) and 
Razoumov & Sommer-Larsen (2010), in order to obtain the relative escape fraction of SNe ionizing photons. This is made possible only by 
using the approximate formula given by equation (13). 
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For each of the stellar spectra we consider, the values we 
find for cr* range from 2.6 x 10"^*^ cm^ to 4.8 x 10"^^ 
cm^, as shown in Table 1. Due to the softer ionizing spec- 
tra of stellar populations at higher metallicity, for higher 
values of the metallicity the average cross section is found 
to be higher as well. Of particular use to us will be the 
average cross section for a stellar metallicity oi Z = 0.05 
Zq , as this is the same population of stars used in the ra- 
diative transfer simulations which we shall draw on in the 
next Section to estimate the escape fractions of SN shock- 
generated ionizing photons (Gnedin et al. 2008; Razoumov 
& Sommer-Larsen 2010; see also Yajima et al. 2011). 

3.2. The Escape Fractions of Stellar and SNe Photons 

In order to compare the escape fractions of stellar and 
SN shock-generated ionizing photons, we first estimate 
their ratio along a given line of sight out of a galaxy into 
the IGM. We then use this result, in conjunction with the 
results of detailed radiative transfer simulations, to esti- 
mate the overall, line of sight-averaged escape fractions of 
ionizing photons. 

3.2.1. Escape Fractions Along a Single Line of Sight 

To begin, we note that the escape fraction of ionizing 
photons along a given line of sight out of a galaxy can be 
approximately expressed a^ 



/esc ^ ^ 



(13) 



where r is the optical depth along the line of sight, which is 
given in terms of the column density Nn of neutral hydro- 
gen and the average cross section a for the absorption of 
ionizing photons by neutral hydrogen and helium. Assum- 
ing that the location of a population of stars in the galaxy 
is the same as that of the SNe that mark their deaths (i.e. 
that the spatial extent of radiating SNe remnants is suf- 
ficiently small) and that the column density iVn along a 
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T ^ 




Nh [cm-2] 

Fig. 1. — Probability distribution of neutral hydrogen column densities A'^h derived from the radiative transfer simulations presented in 
Gnedin et al. (2008) {blue line) and Razoumov & Sommer-Larsen (2010) {red line). The transformation from the escape fraction probability 
distributions presented by these authors is described in Section 3.2.2. Note that, as these authors report, there is a large probability that 
/esc,* = or /esc,* = 1; however, as these values imply A'^h = oo and A'^h = 0, respectively, these points are not shown in this logarithmic plot. 
Nonetheless, they are accounted for in the calculations of the escape fractions of stellar and SN shock-generated ionizing photons presented 
in Section 3.2.2. 



generated photons are larger relative to the time-averaged 
stellar photon escape fractions than we have found here. 

3.2.2. Line of Sight-Averaged Escape Fractions 

To better estimate the relative numbers of hydrogen ion- 
izations caused by SNe and stars, we turn to the results 
of detailed radiative transfer simulations from which the 
escape fraction of stellar photons has been calculated. As 
our result for the ratio of escape fractions is defined only 
along a given line of sight, we must use the results of nu- 
merical simulations which are likewise given along a rep- 
resentative collection of lines of sight. To this end, we 
use the probability distributions for the escape fractions 
of stellar ionizing photons given in figure 2 of Gnedin et 
al. (2008) and figure 9 of Razoumov & Sommer-Larsen 
(2010j3. 

As we would like to generalize our results for the various 
values of cr, given in Table 1, we shall convert the prob- 
ability distributions of escape fractions which are given 
by these authors to the probability distributions of neu- 
tral hydrogen column densities TVh- To do this, we first 
note that from equation (13) we have A^h = -ln(/csc,*)/o'*, 
where /esc,* is the stellar photon escape fraction derived 
from the simulations and cr* is the average cross section 
for photoionization that we have calculated for the same 
stellar spectra used in these simulations, the Z = 0.05 
Zq model from Leitherer et al. (1999). Using this, we 
next convert from the probability distribution dP/df esc,* 
of stellar ionizing photons to the probability distribution 
of neutral hydrogen column densities dP/dNn, as follows: 

^ We note that probability distribution functions similar to those of 
Yajima et al. (2011). 



dP 
dAn 



fci 



dP 



dfc: 



(16) 



The resulting probability distributions dP/dN^i that we 
find for the simulations presented in Gnedin et al. (2008) 
and Razoumov & Sommer-Larsen (2010) are shown in Fig- 
ure 1. These probability distributions contain the basic 
information that we would like to glean from these de- 
tailed hydrodynamics simulations. Using these, we can 
next apply different cross sections cr* and usn to infer the 
probability distributions of the escape fractions of stellar 
and SN shock-generated ionizing photons, and from this 
the overall escape fractions averaged over all lines of sight. 
We note that in our calculation of A^h we have neglected 
the absorption of ionizing photons by dust; however, as 
Gnedin et al. (2008) show, this is a relatively minor effect 
and hence should not impact our results greatly (see also 
Yajima et al. 2011). 

In order to directly compare the probability distribu- 
tions of the two escape fractions /esc,* and /csc,sn, we must 
convert from that for stellar photons dP/df esc,* given in 
Gnedin et al. (2008) and Razoumov & Sommer-Larsen 
(2010) to that for SNe photons. For the stellar spectra 
used in these works, the Z = 0.05 Zq stellar population 
for which we have calculated the average cross section to 
photoionization, this is done as follows: 



dP 



dfc: 



df. 



esc,SN 



csc,SN 



dP 



dfe: 



= 3.7/, 



2.7 

osc,SN 



dP 



dfc: 



(17) 



Razoumov & Sommer-Larsen (2010) have recently been presented in 
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Fig. 2. — Probability distribution functions of escape fractions /esc of ionizing photons from stars {thin lines) and from SNe remnants {thick 
lines), assuming stellar spectra from a population of stars at metallicity Z = 0.05 Zq, formed with an 1 - 100 Mq Salpeter IMF (Leitherer et 
al. 1999). The distribution functions of escape fractions /esc,* of ionizing photons emitted from stars are taken from the results of radiative 
transfer simulations presented in Gnedin et al. (2008) {blue lines) and Razoumov &; Sommer-Larsen (2010) {red lines). The transformation 
of the probability distributions of /esc,* to those of SNe-generated ionizing photons /csc.SN is described in Section 3.2.2. The overall escape 
fractions, averaged over all lines of sight, of SNe and stellar photons are 0.076 and 0.021, respectively, for the results of Gnedin et al. (2008); 
those for the results of Razoumov & Sommer-Larsen (2010) these are 0.80 and 0.22, respectively. The factor of ~ 4 difference between these 
average values of /esc* and /escSN is due to the higher average energy, and so lower cross section for absorption, of the photons generated 
in SNe shocks. Here, this effect is evident in the probability distributions of SNe photons being shifted to higher escape fractions relative to 
those of stellar photons. 




Fig. 3. — The same as Figure 2, but now assuming stellar spectra from a population of stars at metallicity Z = Zq, formed with an 1 - 
100 M0 Salpeter IMF (Leitherer et al. 1999). For this case, the overall escape fractions, averaged over all lines of sight, of SNe and stellar 
photons are 0.076 and 0.011, respectively, for the results of Gnedin et al. (2008); those for the results of Razoumov & Sommer-Larsen (2010) 
these are 0.8 and 0.12, respectively. The factor of ~ 7 difference between the average values for /esc,* and /esc,SN is much larger than in the 
case of the Z = 0.05 Zq spectra shown in Figure 2, due to the lower average energy, and so higher cross section for absorption, of the photons 
from the higher metallicity stars. 
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where we have used equation (15) relating /esc,SN to /esc,*- 
With this conversion from the stehar to the SNe photon 
escape fraction probabihty distribution, we can directly 
compare the two distributions, which we show together in 
Figure 2. As the Figure shows, there is a higher prob- 
ability for larger escape fractions of SN shock-generated 
photons than for stellar photons. In Figure 3, we show the 
probability distributions of /esc,* and /esc.SN for the case 
of a stellar population at solar metallicity, using the ap- 
propriate o", from Table 1. along with the corresponding 
transformation to the probability distribution dP /dfesc,* 
for this average cross section. Comparing the two Figures 
it is evident that the probability distributions for /osc,sn 
and /esc* are much more different in the case of the so- 
lar metallicity stellar population than in the case of the 
lower metallicity stellar population. This is simply due to 
the larger difference between the cross sections it* and ctsn 
for the more metal rich population; because the spectra of 
ionizing radiation from more metal-rich stars is softer, the 
average cross section for photoionization is higher, and so 
the escape fraction of stellar ionizing photons is lower rela- 
tive to the escape fraction of SN shock- generated photons. 

Using our results for the probability distributions of 
both /esc,* and /esc,SN we calculate the overall escape frac- 
tions, averaged over all lines of sight, of SNe and stellar 
photons. For the case of SNe photons, this calculation is 
done as follows: 



< /esc,SN > = 



f 

j osCjSN "7f 
«/esc,SN 



c(fesc,SN 



(18) 



where dP/d/csc,SN is given by equation (17). A similar cal- 
culation is carried out for the line of sight-averaged escape 
fraction of stellar photons < /esc,* >, with the integral 
performed over /esc,* instead of /esc,SN- 

The ratios of the line of sight averaged escape fractions 
that we find for the stellar populations that we consider 
are similar for both of the simulations that we have drawn 
upon, and we present these ratios in the fifth column of Ta- 
ble 1. As shown there, depending on the metallicity of the 
stellar population, we find that the values of /osc.sn/ /csc,* 
vary between ~ 4 and ~ 7. We emphasize that, especially 
for the results of Gnedin et al. (2008), the SNe photon 
escape fraction is significantly lower than what would be 
found by simply applying equation (15) to the overall stel- 
lar escape fraction directly. This is likely due in part to the 
fact that the simulation results show relatively large prob- 
abilities for escape fractions of either /esc,* = 1 or /esc,* = 
0; because these limiting values will yield the same escape 
fractions for SN generated photons, via equation (15), the 
overall line of sight-averaged escape fractions of stellar and 
SNe photons are not as different as they are along individ- 
ual lines of sight for which < /esc,* < 1- 

Finally, we note that our results are senstively depen- 
dent on the probability distribution functions that we have 
adopted from Gnedin et al. (2008) and Razoumov & 
Sommer-Larsen (2010). While the detailed simulations 
which these authors carry out in attaining these results 
are some of the most sophisticated and complete to date 
(including radiative feedback from young stars and both 
mechanical and chemical feedback from SNe), the resolu- 
tion of the simulations is limited and the structure of the 
interstellar medium at sub-resolution scales will have an 



impact on the actual escape fractions of ionizing photons. 
Despite this uncertainty which plagues all escape fraction 
calculations, the results presented by these authors allow 
us to make a first estimate of the relative role of SNe in 
cosmic rcionization. Future simulations at higher resolu- 
tion will ideally include the ionizing radiation from both 
stars and SNe explicitly, thereby allowing for improved es- 
timates of the relative role of each of these sources. 

4. THE RELATIVE CONTRIBUTIONS OF STARS AND SNE 
TO REIONIZATION 

Working from our results from Sections 2 and 3, we may 
now compare the relative contributions of stars and SNe to 
reionization. The ratio of the number of hydrogen ioniza- 
tions in the IGM due to SNe to that due to stars is given by 
< /oscSN >Qosc,sn/< ./csc,* >Qcsc,*. Usiug cquatious (9) 
and (10); along with the values of and iVgN from Lei- 
therer et al. (1999) and the ratios of line of sight-averaged 
escape fractions, each given in Table 1, we estimate this 
ratio as 



< /esc,SN > QsN ^ Q ]^ I 



^ /esc,* ^ 



Vl053g 

X 

/ < /esc,SN > \ / < /esc,* > \ gx 
V 0.4 y ^ 0.1 / 

The escape fractions appearing here are the overall line 
of sight-averaged values, normalized to typical values. In 
particular, we have normalized the ratio of escape fractions 
to < /esc,SN >/< /esc,* > = 4, roughly the ratio found in 
the last Section for a stellar population at sub-solar metal- 
licity. 

The right-most column of Table 1 shows the values of 

this ratio for all of the stellar populations we consider. The 
values that we find for the contribution of SNe to reioniza- 
tion, relative to the contribution from stars, ranges from 
~ 0.05 to ~ 0.12, with higher relative contributions for the 
more metal-rich stellar populations. As discussed in Sec- 
tion 3, this is due largely to the higher ratios of escape frac- 
tions < /esc,SN >/< /csc* > for the more metal-enriched 
populations, owing to their softer stellar spectra. However, 
the different rates of SN production A^sn and stellar ioniz- 
ing photon generation also have an effect. Importantly, 
we note that the relative contributions to rcionization scale 
linearly with the energy per supernova E'sn ; while we have 
normalized our results to £Jsn = 10^^ erg, if SNe are more 
energetic than this in the early universe, then the relative 
contribution of SNe to reionization may be considerably 
higher. 

Finally, we note that our calculations of the quantities 
appearing in Table 1 have been done drawing on numerous 

approximations. We have nonetheless chosen to report the 
values of these quantities, in cases, to two or three deci- 
mal places, in order to highlight the differences between 
the stellar models at different metallicities that we have 
considered. 

5. THE CONTRIBUTION OF SUPERNOVAE-GENERATED 
COSMIC RAYS TO REIONIZATION 
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Table 1 

Relative contributions of stars and Supernovae to reionization 



Mctallicitv'^ 
[Zq\ 


[105?s-i] 


iVsN 

[10-2 M-1] 


[10-18 cm2] 


Ratio of escape fraetions"^ : 

< /esc.SN >/ < /esc,* > 


Relative eontribution of SNe: 

< /esc.SN >(3sn/ < /esc,* >Q*'^ 


1 


2.2 


1.91 


4.8 


6.8 


0.12 


0.4 


2.5 


1.86 


2.9 


4.2 


0.062 


0.2 


2.8 


1.91 


2.8 


4.1 


0.056 


0.05 


3.2 


1.95 


2.6 


4.0 


0.048 



"The metallicity of the stellar population. Here we take the solar iiietallicity to be Zq = 0.02. 
''The number of ionizing photons per second produced for a constant star formation rate of 1 My yr^^. 
^Estimated from the column density probability distributions derived from the results presented in Gnedin et al 
Sommer-Larsen (2010), as described in Section 3.2.2. 

'^Note that the relative contribution of SNe scales linearly with the average SN energy Eg^ ; here we assume that EgN 



(2008) and Razoumov 
lO^i erg. 



Note. — All stellar populations modeled here are assumed to have a Salpeter IMF, with stars formed in the mass range 1 - 100 Mq. The 
values of Q*, ATgN, and cr« were computed from the models of these stellar populations presented in the Starburst99 database (Leitherer et al. 
1999), assuming a constant star formation rate. Finally, we have assumed a constant cgN = 0.7 x 10-^* cm^, independent of metallicity, as 
discussed in Section 3.1. 



An additional process driven by SNe which leads to ion- 
ization of hydrogen is the production of cosmic rays, and 
for completeness here we consider their contribution to 
reionization. We may find an upper limit to the total 
number of ionizations that cosmic rays can cause in the 
IGM, by making the simple assumption that all the en- 
ergy injected into cosmic rays by SNe is deposited into the 
IGM via ionizations of hydrogen. Following the discussion 
of cosmic ray production in the early universe given by 
Stacy & Bromm (2007), we shall assume that 10 percent 
of the energy i?sN of SNe goes into cosmic rays (e.g. Ru- 
derman 1974); then, noting that each ionization removes 
50 eV of energy from the impinging cosmic ray (Spitzer & 
Tomasko 1968), we find the maximum number of ioniza- 
tions that can be produced per SN to be ~ 10^°. Relating 
this to the star formation rate, in order to compare the 
maximum rate (assuming an escape fraction of cosmic rays 
of /esc.CR = 1) of ionizations Qcr by cosmic rays to the 
rate of ionizations due to stellar and SN shock-generated 
photons, we have 



Qcr ^ 4 X 10 



50 



M© yr- 



10-2 M^^ j \ 1051 erg 



(20) 



We emphasize that this is an extremely hard upper limit, 
as in general the IGM is optically thin to cosmic rays, es- 
pecially to those with higher energies. To illustrate this, 
we can make a simple conversion of the result presented 
in figure 1 of Stacy & Bromm (2007) to find the penetra- 
tion depth dcK of a cosmic ray with energy ecR traveling 
through a completely neutral general IGM at a redshift z: 



dcK ~ 10 



l + z 
10 



(l^) ^P^' 



where we have applied a simple power-law fit appropri- 
ate for low cosmic ray energies (i.e. ecR. ^ 10^ cV), and 
the distance is in physical units. Comparing this to the 
physical distance that a cosmic ray moving at close to the 
speed of light would travel from its generation at 2; < 20 in 
a SN explosion to the end of reionization at 2; ~ 6, which 
is < 200 Mpc, we see that only cosmic rays with ecR ^ 
lO*" eV will have time to deposit a large portion of their 
energy into ionizing the ICM. Although the energy distri- 
bution of cosmic rays in the early universe is not known, 
it is likely that the bulk of the total energy in cosmic rays 
generated before reionization is carried by cosmic rays at 
higher energies that this (e.g. Jasche et al. 2007; Stacy & 
Bromm 2007). Therefore, comparing Qck to Q* and Qsn, 
we conclude that the contribution to reionization from cos- 
mic rays generated in SNe is likely to be well below that 
from either stellar sources or from the radiation generated 
in SNe shocks. 

6. IMPLICATIONS FOR REIONIZATION BY GALAXIES 

In the previous Sections we have estimated the relative 
contributions of SNe and stars to reionization. We now 
consider the implications of our results for the population 
of galaxies responsible for reionization. In particular, we 
highlight how constraints on the luminosity function of 
reionizing galaxies are strengthened by accounting for the 
contribution of SNe. 

Following Wilkins et al. (2011), we assume that the 
number of ionizing photons produced by a galaxy of lumi- 
nosity L1600 at 1600 A is = Dion x IQi^Lieoo, where the 
units of iieoo is erg Kr^. To facilitate a comparison 
with this work, we use for the fraction of ionizing photons 
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Fig. 4. — The limiting magnitude down to which the luminosity function needs to be integrated to achieve reionization at ^ = 7 for a given 
ratio of clumping factor C to escape fraction of stellar photons /esc,*- Solid lines show results for stellar photons only assuming the fits to 
the luminosity function presented in Wilkins et al. (2011). Faint-end slopes a are increasingly steep, going from left to right. Short-dashed, 
dot-dashed and long-dashed lines show the corresponding cases including the contribution from SN shock-generated photons to reionization. 
Note that this Figure differs from the similar figure 17 of Wilkins et al. (2011), who used Dion = 1-47 instead of Dion = 1.2 (Wilkins, private 
communication) . 



-Dion = 1-2. This is appropriate for a stellar population 
with metallicity Z = Zq, with stars having formed with a 
Salpeter IMF at a constant rate for 100 Myr (Leitherer et 
al. 1999), as these authors consider. The rate of ionizing 
photon production per unit comoving volume that is re- 
quired to maintain reionization is given by (Madau et al. 
1999; Wilkins 2011): 

Qreq = 2.51 X 10^^C(1 + zf S'^ Mpc'^ . (22) 

Here C is the clumping factor of the IGM. Comparing the 

above expression for q^^q to the rate of photon production 
by both stars and SNe, taking into account our result from 
Table 1 that 0.12/esc,*Q* — /csc,snQsn, we integrate over 
the galaxy luminosity function to find the following for the 
ratio of clumping factor to escape fraction required for the 
universe to remain reionized: 

= 1.12 X ^O'^^T^^i^ [ Lieoo(t){Lieoo)dLieoa ■ 

/esc,* ~^ ^) J 

(23) 

For the luminosity function (j) we use the z = 7 best- 
fit values reported in Wilkins et al. (2011). Given 
the uncertainty in the faint-end slope of the luminos- 
ity function (e.g. Ryan et al. 2007; Khochfar et al. 
2007), we predict the ratio of C//* for three different 
choices of faint-end slopes {a = —1.5, —1.7, —1.9; ^* = 
0.00126, 0.00106, 0.00072; Mi*6po = -19.8,-19.9,-20.1). 
Typical values for the clumping factor and escape frac- 
tion are C ^ 5 (e.g. Pawlik et al. 2009) and /esc,* ~ 0.2 
(e.g. Ciardi & Ferrara 2005), respectively, which yields 
C//esc,* ~ 25. As shown in Figure 4, such high values 



are only reached if the luminosity function of galaxies is 
integrated down to a limiting magnitude of -Mieoo ~ —8 
and if the faint-end slope of the luminosity function is as 
steep as a = —1.9. In the case of only stellar photons the 
limiting magnitude is fainter by several magnitudes, show- 
ing that the contribution to reionization from stars in low 
luminous galaxies is comparable to the one from SNe in 
massive galaxies. 

Finally, we note that recent work suggests that the star 
formation rates of high rcdshift galaxies are not constant, 
as assumed in the calculation of the ionization parame- 
ter Dion used above, but instead are increasing with time 
(e.g. Finlator et al. 2011; Khochfar & Silk 2011). In 
this event, I?ion should be considerably higher (Wilkins et 
al. 2011), and thus the contribution from low luminosity 
galaxies may be even lower than we find here. 

7. DISCUSSION AND CONCLUSIONS 

We have used state of the art models of stellar popu- 
lations and high velocity shocks, along with the results 
of detailed radiative transfer simulations, to evaluate the 
contribution of SNe to the reionization of the universe. We 
have found that, although the number of ionizing photons 
emitted from SN shocks is well below that emitted from 
stars, due to the harder spectrum of the ionizing radia- 
tion from SN shocks the fraction of such photons which 
escape galaxies is in general much larger than the escape 
fraction of stellar ionizing photons. Related to this, ow- 
ing principally to the softer spectrum of stars of higher 
metallicity, we find that the relative contribution of SNe 
to reionization is metallicity dependent; in particular, for 
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stellar metallicities between 0.05 Zq ^ Z < Zq, SNe are 
responsible for between ~ 5 and ~ 12 percent of the num- 
ber of ionizations caused by stellar sources. Also, while 
these numbers are for line of sight-averaged escape frac- 
tions, we have furthermore found that, along a given line 
of sight out of a galaxy from a star-forming region, for 
hydrogen column densities iVn ^ 10^^ cm~^ the contri- 
bution to reionization from SNe is larger than that from 
stars. Finally, we have shown that the contribution from 
SNe allows reionization to be completed by galaxies a few 
magnitudes brighter than in the case in which only stellar 
sources contribute. 

We emphasize that we have found a conservative esti- 
mate of the contribution from SNe, for the following rea- 
sons. Firstly, we have included only the ionizing radia- 
tion emitted from SN shocks in the Sedov- Taylor phase; 
while we find the radiated energy to already be ~ 40 per- 
cent of the total SN energy at this stage, there may be 
an additional contribution of ionizing photons from SNe 
in the later, radiative phase if the upstream medium is 
already ionized by other sources (ShuU & Silk 1979; see 
also Section 2.1.1 here). Secondly, in our modeling we 
have not accounted for the fact that the majority of SN 
shock-generated photons are emitted after the stellar ion- 
izing radiation from a given population of stars; as the 
radiation from the stars acts to ionize and rarify the gas 
along lines of sight out of a galaxy, when SNe explode the 
column density of neutral hydrogen A^h will likely be lower 
than the average column density seen by stellar photons. 
We also note that, as we find that the number of ion- 
izations caused by SN scales linearly with SN energy, if 
this is higher or lower than our fiducial value of -Esn — 
\q5i g]-g during the epoch of reionization, then SNe may 
have contributed more or less, respectively, to reionization 
than we have found here. Related to this, if the stellar 
IMF is more top-heavy during reionization than it is in 
the Galaxy today, as some recent studies would suggest 
(e.g. van Dokkum 2008; Dave 2008; Gunawardhana et al. 



2011), then this may also lead to an enhancement in the 
contribution from SNe, as in this case a higher fraction of 
stars may explode as SNe and these SNe may also be more 
energetic. 

Beyond contributing to the ionization of hydrogen in the 
IGM, the enhanced photoheating rate due to the contri- 
bution of ionizing photons from SNe likely has additional 
effects. For instance, a higher heating rate of the gas con- 
tributes to the suppression of low-mass galaxy formation 
in reionized regions of the universe (e.g. Dijkstra et al. 
2004; Hambrick et al. 2010). It is also likely to lead to a 
reduction in the clumping factor C of the IGM, which has 
bearing on the rate at which reionization occurred; how- 
ever, we note that for heating rates comparable to what 
would be expected for SN-generated ionizations Pawlik et 
al. (2009) find only a small decrease in the clumping fac- 
tor. 

While we have found SNe to be of secondary impor- 
tance compared to stars as contributors to reionization, 
they undeniably provide an enhancement in the rate at 
which reionization occurs above what is found in models 
accounting for only stellar sources (e.g. Sokasian et al. 
2003; Iliev et al. 2006; Choudhury & Ferrara 2007; Trac 
& Cen 2007; Raicevic et al. 2011). In this, accounting for 
the effect of SNe can aid in explaining how the universe was 
reionized by 2; ^ 6. Using the results we have presented in 
Table 1, it is straightforward to include the contribution of 
SNe in reionization calculations that already include the 
contribution from stars. 
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APPENDIX 

AN APPROXIMATE EXPRESSION FOR THE ESCAPE FRACTION ALONG A SINGLE LINE OF SIGHT 

Here we justify the expression that we have used for the fraction of ionizing photons which escape along a single line 
of sight, equation (13). The exact expression for the escape fraction of ionizing photons along a given line of sight with 
column density of hydrogen A^h is given by 

r J13.6eV -fc-, ' / A 1 \ 

/<^^c,* - — F(E^) ^j^ ' ^^^) 

Jl3.6eV "-^7 

where for concreteness we have chosen to treat the case of stellar photons; a similar treatment to that given here can 
easily be given for the case of SNe photons. In order to show that equation (13) is approximately equal to this expression, 
we expand the exponential in the numerator, which yields 



/i?.6eV^[l-^Ha.(^;,) ] dE. 

rex. F{E^) ,p, 
Jl3.6eV E-, "-^7 



Slightly rewriting this, we have 
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(A2) 



/e-* ^ 1 - ■-■7; . (A3) 

Jl3.6eV E^ ""^l 

Noting that the cross section here is that due to absorption of ionizing photons by both hydrgen and helium, as discussed 
in Section 3.1, we have (T*(£'-y) = amiE-^) + cFYif.i{E^). Then, with equation (12), we obtain 

/esc* - 1 - ^H^* ~ e-^"'^' , (A4) 

which is eqtiation (13) for the case of stellar photons. This demonstrates that equation (13) is valid as an approximate 
expression for the escape fraction. 



